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Helicenes are inherently chiral polyaromatic molecules composed of all-ortho fused benzene rings 
possessing a spring-like structure. Here, using a combination of density functional theory and tight- 
binding calculations, it is demonstrated that controlling the length of the helicene molecule by 
mechanically stretching or compressing the molecular junction can dramatically change the elec¬ 
tronic properties of the helicene, leading to a tunable switching behavior of the conductance and 
thermopower of the junction with on/off ratios of several orders of magnitude. Furthermore, control 
over the helicene length and number of rings is shown to lead to more than an order of magnitude 
increase in the thermopower and thermoelectric figure-of-merit over typical molecular junctions, 
presenting new possibilities of making efficient thermoelectric molecular devices. The physical ori¬ 
gin of the strong dependence of the transport properties of the junction is investigated, and found 
to be related to a shift in the position of the molecular orbitals. 


INTRODUCTION 


Interest in single-molecule junctions arose from 
their originally proposed role as functional ele¬ 
ments in electronic devices [1]. In recent years 
it has been realized that molecular junctions 
(MJs) can potentially fulfill additional function¬ 
alities ranging from opto-clectronics and spin- 
tronics to thermoelectric energy conversion [2- 
4], which is the focus of this work. The ther¬ 
mopower and the thermoelectric figure of merit 
(FOM) are measures of a junction’s ability to 
convert a temperature difference into electric 
power, and MJs were suggested as candidates 
for efficient and high-power thermoelectric de¬ 
vices [5-14] due to their low dimensionality, 
large versatility and low thermal conductance 
[4, 15, 16]. While there have already been several 
demonstrations of thermoelectric energy conver¬ 
sion in single-molecule junctions [17-24], unfor¬ 
tunately the typical thermopower is rather small 
(S' ~ 5 — 50 /jlV/ K), much smaller than that 
of commercially available semi-conductor-based 
thermoelectrics. Increasing the thermopower of 
molecular junctions is an essential step if MJs 
are ever to become relevant energy conversion 
technologies. The natural question that arises 


in this context is then: is there a way to in-situ 
increase the thermopower, S, and the thermo¬ 
electric FOM, ZT, of MJs? 

In this letter, we propose to focus on hclicene- 
based molecular junctions (HMJs), which can 
be used as a tool to help answer the above 
question. [n]Helicenes are aromatic com¬ 
pounds in which n benzene rings or other (het- 
ero)aromatics are angularly annulated to give 
helically-shaped molecules with a spring-like 
structure (Fig. 1(a)). The chemistry of he¬ 
licenes has been systematically explored for more 
than half a century (for recent reviews see, 
e.g. [25-28]), and a renewed and growing inter¬ 
est in helical aromatics is clearly visible within 
the last decade. Although the number of ap¬ 
plications of helicenes is so far rather limited, 
they have already found astonishing applications 
to enantioselective organo- or transition metal 
catalysis [29, 30], molecular recognition, [31, 32] 
self-assembly, [33] nonlinear optical materials, 
[34-36] chiral materials, [37-39] and chiroptical 
electronic devices. [40-42] Furthermore, as a 
first step towards realizing single-molecule junc¬ 
tions, helicenes have been placed on surfaces and 
probed with atomic force and scanning tunneling 
microscopes. [43-50] 



2 



2,15-diaza[6]helicene (1) 2,17-diaza[7]helicene (2) 2,19-diaza[8]helicene (3) 



2,27-diaza[12]helicene (7) 2,29-diaza[13]helicene (8) 2,31-diaza[14]helicene (9) 



FIG. 1. Schematic illustration of the system studied: (a) 
Chemical structure of diaza[n]helicenes (l)-(9) used in the 
DFT and tight-binding calculations. Chirality issues were not 
a subject of this study, however, all the investigated systems 
were of (M) helicity. (b) 2,31-Diaza[14]helicene (9) molecule 
connected to two gold(l,l,l) contacts, r denoting the dis¬ 
tance between the electrodes, (c) Graphic representation of 
the tight-binding model of the helicene-based molecular junc¬ 
tion. The tight-binding parameters for nearest neighbor in¬ 
teraction ( t ), inter-stack coupling (fi) and molecule-electrode 
level broadening V are noted. The four points in (c) mark 
the contact points of the rings with the electrode used for the 
tight-binding model. 


We suggest that the distance r between the two 
metallic electrodes of a hclicene molecular junc¬ 
tion (Fig. 1(b)) can be used as an experimen¬ 
tal control parameter that can be tuned to al¬ 
ter and probe the properties of the HMJ. Due to 
the spring-like structure of the hclicene molecule, 
when the junction is compressed, i.e., the dis¬ 
tance between the electrodes of the junction is 
reduced, the couplings between the carbon or¬ 
bitals on the benzene rings increases, opening 
additional transport paths for the electrons and 
coherently changing the electronic properties of 
the MJ. Regulation of the electrode distance in 
molecular junctions has already been introduced 
as a means to control and measure the molec¬ 
ular forces and couplings [51-54] (interestingly, 
fullcrene C60 was already demonstrated to mini¬ 
mize its resistance and become almost transpar¬ 


ent for tunneling electrons if the molecule is suf¬ 
ficiently squeezed by an STM tip). [55] How¬ 
ever, in the studies mentioned above, the domi¬ 
nant origin for changes in transport is the change 
in the molecule-electrode coupling, while the 
molecular electronic structure remains largely 
unchanged. In a nonplanar molecule, on the 
other hand, pulling or stretching should affect 
both the molecular orbitals and the couplings. 
Nevertheless, there is a lack of suitable models 
of nonplanar molecules to systematically investi¬ 
gate the stretch- or compression-dependent sin¬ 
gle molecule conductance. 

We demonstrate our conceptual experiment 
by using density-functional theory (DFT)-based 
transport calculations to extract the transmis¬ 
sion function of a HMJ for different values of 
the inter-electrode distance, r. To demonstrate 
how the hclicene electronic structure depends on 
r, we then use a tight-binding (TB) model to fit 
the DFT data. Next, the TB model and param¬ 
eters are used to predict the thermopower and 
ZT of hclicene molecular junctions as a func¬ 
tion of r and helicene length, demonstrating that 
these can be optimized to obtain maximal ther¬ 
moelectric efficiency. It is important to point 
out early on that we are aiming to generate 
qualitative rather than quantitative predictions, 
and therefore the use of a TB approximation is 
justified. However, even when using the first- 
principle DFT calculation, the predictions can 
only be qualitative due to the lack of appropri¬ 
ate functionals and benchmark DFT calculations 
for helicenes. 


METHODS 

Molecular system under study 


We studies diazahelicenes with the hetero¬ 
atoms in position 2 (symmetrically on both 
sides). These systems were selected for three rea¬ 
sons: (i) they are fully aromatic and highly con¬ 
ductive, (ii) they are synthetically accessible and 
were previously prepared by some of the authors 
[56] and (iii) they provide for reasonably good 
binding to gold surfaces via the nitrogen func¬ 
tionalities. [57-59] For this reason, azamolecules 
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such as 4,4’-bipyridine are often used in single 
molecule conductance experiments, for example 
break junction measurements. [60] Based on our 
DFT calculations, position 2 in helicenes appears 
more suitable for binding than, for example, po¬ 
sition 1, for steric reasons. 

To achieve conductance and thermopower tun¬ 
ing and switching we used two additional pa¬ 
rameters. One is helicene length ([n]) which 
was varied between 6 and 14 aromatic rings and 
the other is helicene stretch/compression (Ar). 
Out of the 9 total molecular systems of differ¬ 
ent number of rings, 2,15-cliaza[6]helicene (1) 
to 2,31-diaza[14]helicene (9), only one was com¬ 
pletely investigated by DFT methods. This was 
the 2,21-diaza[9]helicene (4) (Fig. 1(a)) from the 
middle of the series. Here we calculated all the 
stretched and compressed geometries, energies, 
and transmission spectra. Due to the system 
size, the calculations were very lengthy and cal¬ 
culation of all 9 molecular structures would be 
impractical. This was the motivation to use a 
simple tight-binding model to study the remain¬ 
ing systems at much lower computational cost. 
The tight-binding model was, however, found to 
reproduce the DFT data surprisingly well. Chi¬ 
rality issues were not a subject of this study, all 
the investigated systems were of (M) hclicity. 


Density-functional-theory calculation 

We considered a 2,21-diaza[9]-helicene(4) 
placed between two Au(lll) electrodes. Prior 
to any periodic DFT calculations, all molecu¬ 
lar structures were optimized in vacuum at the 
B97D/cc-pVDZ level of theory using the G09 
package. [61] QuantumWise Atomistix toolkit 
(ATK) [62-64] was then employed in all subse¬ 
quent calculations described below. First, the 
HMJ device was constructed as a ” two-probe 
system” within the Virtual Nano Lab (VNL) 
module of ATK. Here, the HMJ consists of 
three regions: The left electrode (bulk Au), the 
so called central region and the right electrode 
(bulk Au), cf. Fig. 2 top panel. In the central 
region the helicene molecule is sandwiched be¬ 
tween six surface layers of Au(lll), i.e. with 


three Au layers at both sides of the molecule. 
The unit cell size of the left and right elec¬ 
trodes in the plane parallel to the gold surface 
was 5x5 Au atoms. The entire two-probe sys¬ 
tem was then optimized and the electrode sep¬ 
aration relaxed. This resulted in our reference 
equilibrium (Ar = 0, corresponding to an elec¬ 
trode separation of ~ 8.5 rA) geometry from 
which we started the stretching and squeezing 
by just adjusting the electrode separation in 0.5 
rA steps. The total compression attempted was 
3.5 rA while the stretching continued until the 
molecule broke off the electrodes at 7.0 rA elon¬ 
gation. Only in the optimal configuration the 
unit cell was allowed to fully relax, so only lim¬ 
ited deformations of the gold electrodes were al¬ 
lowed. The system was re-optimized at each step 
with the electrode separation constrained. DFT- 
PBE/SZP level of theory with 2 x 2 x 50 k-point 
sampling was used for all optimizations while 
DFT-PBE/DZP with 3x3x100 k-point sampling 
was used for transmission spectra evaluation. 
To check the numerical accuracy, optimizations 
with larger k-point sampling (3 x 3 x 100) were 
also tried with no difference in the resulting ge¬ 
ometries. The bottom panel of Fig. 2 shows the 
total energy E as a function of Ar. The forces 
required for stretching and squeezing can be es¬ 
timated from the energy curve to be ~ 1 nN and 
5 nN, respectively. 


Transport calculation within the tight-binding model 

The DFT results indicate that helicenes tend 
to lay flat on the electrode surface whenever pos¬ 
sible [26, 65]. At least one aromatic ring (but 
more than one for most geometries) tend to form 
contact with the electrodes for the electrode sep¬ 
arations used in the TB and thermopower calcu¬ 
lations, cf. Fig. 2. Of course, the contact weak¬ 
ens as we stretch the molecule out but a four- 
point contact model appeared to be applicable 
throughout the stretched/compressed structures 
used in this paper (Fig. 1(c)). 

Once the Hamiltonian and self-energies are de¬ 
termined (Eq. 1-2), all the transport proper¬ 
ties can be determined within the Landauer for- 
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FIG. 2. Top panel: Two-probe model of the 2,21- 
diaza[9]helicene (4) HMJ used in the DFT geometry optimiza¬ 
tions and transport calculations. Compressed (top), optimal 
(center) and stretched (bottom) structures. The green boxes 
in the optimal geometry indicate the left and right electrode 
regions, while the red box encircles the central region. Bottom 
panel: total energy E as a function of A r. 


malism [66-68]. The Green’s functions are de¬ 
termined via G r,a = (E — T-L m + F r,a )~ 1 , where 
E r,a = -f Eg“, corresponding to the self¬ 
energies of the top and bottom electrodes. 
The transmission function is given by T(E) = 
Tr (E r T G r 'E c j 3 G a ), and the transport coefficients - 
the conductance G, the thermopower S, and the 
thermal conductance /t - are determined within 
the Landauer formalism as G = e 2 Lo, S = 


Li/(eTL 0 , k = [l 2 — /T, where T is the 

temperature (room temperature is considered 
hereafter) and L n = — 4 J dET(E)(E — p) n f^ 
are the Landauer integrals, with h being Planck’s 
constant, /i the chemical potential of the elec¬ 
trodes (set as the zero of energy hereafter), and 
f(E) the Fermi-Dirac distributions. The ther¬ 
moelectric figure of merit ZT is given by ZT = 
GS 2 
k/T ' 

RESULTS: TUNING CONDUCTANCE AND 
THERMOPOWER 

The starting point is thus a DFT-based cal¬ 
culation of the transmission through a he- 
licene molecular junction, formed by 2,21- 
diaza[9]hclicene (4) and two gold(l,l,l) elec¬ 
trodes. Chirality issues were not a subject of 
this study, however, all the investigated systems 
were of (M) helicity. The transmission curves 
were calculated with the combined DFT-non- 
equilibrium Green’s function approach [62-64], 
Fig. 3(a) shows the transmission curves T(E) 
as a function of energy E for different inter- 
electrode distances r, incrementally changed 
about A r = —3, —2.5, —2,..., 1,1.5rA , mea¬ 
sured from the relaxed inter-electrode distance 
r (it is important to point here that Ar = 0 
is the distance where the DFT calculation finds 
the energy minimum, but may not correspond to 
any experimental distance, since in any experi¬ 
mental setting the distance will be controlled by 
an external force). The most striking feature of 
Fig. 3(a) is the apparent closing of the transmis¬ 
sion gap between the highest occupied and low¬ 
est unoccupied molecular orbitals (HOMO and 
LUMO, respectively) upon compression of the 
molecular junction, which transforms the MJs 
electronic nature from ’’insulating” to ’’metal¬ 
lic” , by which we mean resonant tunneling (close 
to a molecular level). Consequently, the conduc¬ 
tance changes by more than three orders of mag¬ 
nitude. In Fig. 3(b) the conductance (in units 
of the conductance quantum Go = 2 e 2 /h) and 
the thermopower are plotted as a function of the 
distance r. The conductance shows a switching 
behavior with an ”on”/’’off’ ratio of more than 
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three orders of magnitude (inset of Fig.3(b)), 
suggesting that HMJs can perform as reversible, 
nanoscale molecular switches. [54, 69-77] The 
thermopower shows a sign-change, a direct ex¬ 
perimental prediction, this is an important fea¬ 
ture, since a sign-change of the thermopower im¬ 
plies a change in the transport mechanism of the 
junction [4], yet such a change is hard to demon¬ 
strate in situ. 


(ii) the TB model allows for interpretation of 
the DFT calculation in terms of a simple physi¬ 
cal picture, which is in many cases hard to build 
from the output of the DFT calculation (typ¬ 
ically energies, orbitals and transmission), (iii) 
additional effects (such as dephasing, see below) 
cannot be introduced within the DFT calcula¬ 
tion. The hclicene molecule TB Hamiltonian is 
given by 



Ar [A] 


FIG. 3. Transmission and thermopower properties of HMJ 
formed by 2,21-diaza[9]helicene (4) and gold(l,l,l) elec¬ 
trodes:^) Transmission T(E) as a function of energy E for dif¬ 
ferent inter-electrode distances incrementally changed about 
Ar = — 3, 2.5,-2,..., 1, 2rA (measured from the relaxed inter¬ 
electrode distance), (b) Conductance (blue circles) and ther¬ 
mopower (yellow diamonds) as a function of Ar. Inset: Con¬ 
ductance on a logarithmic scale, demonstrating the order of 
magnitude change in conductance with changing Ar. 


To understand the origin of this change in the 
hclicene electronic structure, we construct a TB 
description of the system. The use of TB models 
to complement DFT calculations is rather com¬ 
mon (see, e.g. [78-81]), and is useful for several 
reasons: (i) DFT calculations may become ex¬ 
tremely computationally expensive and lengthy, 


U M = X/ol^l I*)01 > 

* <V'> { i,j} 

( 1 ) 

where \i) represents the orbital located on the i- 
tli atom in the hclicene, eo is the orbital energy, 
t is the nearest-neighbor hopping matrix ele¬ 
ments (with (i, i') representing nearest-neighbor 
orbitals), and t\ represents hopping matrix el¬ 
ements between sites that reside vertically one 
above the other on adjacent turns of the helix, 
so-called inter-stack hopping (with {i,j} repre¬ 
senting vertically-neighboring orbitals), t and t\ 
are graphically depicted in Fig. 1(c). The elec¬ 
trodes are treated within a wide-band approxi¬ 
mation [82] through a self-energy term, 

= |*t,b)(*t,b| , (2) 

*T,B 

where r(a) stands for retarded (advanced) self 
energy, T and B stand for top- and bottom- 
electrodes, respectively, and V TB are the cor¬ 
responding electrode-induced level-broadenings 
(we set T^ = T# = T for simplicity hereafter). 
The indices It,b represent the sites on the he- 
licene molecule that are in contact with the elec¬ 
trode (see Fig. 1(c)). We assume the edge ben¬ 
zene rings of the hclicene couple with the elec¬ 
trodes via the four external orbitals (that is, the 
edge rings lie almost flat on the electrodes). We 
found that changing the contact configuration 
does not qualitatively change our results. 

To obtain the TB parameters, parameters 
£o,t,ti and T are tuned for a best fit between 
the transmission functions as obtained from the 
tight-binding and from the DFT for a HMJ 
(formed by 2,21-diaza[9]hclicene (4) and gold 
(1,1,1) electrodes) at Ar = 0. In the top in- 
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set of Fig. 4 the transmission curves of the DFT 
(solid blue line) and TB (dashed red line) cal¬ 
culations are shown, with the best fit yielding 
the tight-binding parameters e 0 = —0.61 eV, 
t = 1.98 eV, t\ = 0.279 eV and F = 1.78 eV. 
The TB curve fits the DFT data very well, even 
with such a small number of parameters and 
especially around the HOMO and LUMO reso¬ 
nances and below the HOMO level (note that for 
transport calculations, the important fit region 
is within fcgT from the Fermi level). 

To capture the HMJ electronic behavior un¬ 
der compression or stretching, it is necessary to 
know the r-dependence of the TB parameters. It 
is safe to assume that the orbital energy e 0 and 
the nearest-neighbor hopping element t do not 
change substantially when the helicene molecule 
is stretched, and that most of the change occurs 
in t\ and T. Based on knowledge that the depen¬ 
dence of t\ and r on r originates from changes 
in the overlap of wave-functions, a reasonable as¬ 
sumption is that they depend exponentially on r, 

ti(r) = h exp , T(r) = F exp . The 

values of £ t and £p are found by fitting the TB 
transmission at the Fermi energy T(ep) (which 
at low temperatures is proportional to the con¬ 
ductance [67, 68]) to the DFT data. In Fig. 4, 
the conductance is plotted as a function of the 
the change in the inter-electrode distance, Ar, 
obtained from the DFT data (filled gray squares) 
and the TB calculation (solid blue line). This 
fit yields the values = 1.33rA and £p = 1-2 
rA . To demonstrate that it is necessary to in¬ 
clude the change in ti, the dashed orange line 
of Fig. 4 shows T(e^) as a function of r for the 
case where only F is dependent on r and t\ is 
independent of r (i.e. where the transport de¬ 
pendence on r comes only from the change in 
the contacts, while the electronic structure of 
the molecule itself is unchanged). The striking 
difference between the dashed orange line and 
the DFT calculation (filled gray squares) clearly 
shows that the internal electronic structure of 
the molecule is changed upon stretching (via a 
change in t\). The bottom inset of Fig. 4 shows 
the thermopower S' as a function of r based on 
the DFT data and TB calculation (solid blue 
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FIG. 4. 2,21-Diaza[9]helicene (4): Conductance at the Fermi 
energy level as a function of the change A r in the inter¬ 
electrode distance measured from the relaxed configuration, 
based on the DFT data (filled gray squares) and on the TB 
calculation, assuming that both t\(r) and T(r) depend expo¬ 
nentially on r (solid blue line). The dashed orange line shows 
the same for a TB calculation where ti is independent of r 
(see text for parameter values). Top inset: the transmission 
as a function of energy for the DFT data (solid blue line) and 
the TB calculation (dashed red line). Bottom inset: Ther¬ 
mopower S as a function of Ar based on the DFT data (filled 
gray squares) and TB calculation (solid blue line). 

line). While some shift is visible, the overall or¬ 
der of magnitude and trend are similar, which is 
impressive considering the small number of pa¬ 
rameters in the model and the fact that ther¬ 
mopower is extremely sensitive to small varia¬ 
tions in the parameters of the molecular junction 
[78, 83] (note that the thermopower is approx¬ 
imately proportional to the logarithmic deriva¬ 
tive of the transmission function). 

To demonstrate the potential of HMJs as 
switches, in Fig. 5 the I — V curve of a [9]he¬ 
licene HMJ is plotted for two values of stretch¬ 
ing distance, r = —2rA (solid blue line) and 
r = 2rA (dashed yellow line). The prominent 
dissimilarity between them shows that the HMJ 
can be mechanically switched from a ’’metallic” 
(”on” state) to an insulating (’’off”) state. The 
bottom-right inset shows the same on a log scale, 
to more clearly show the four orders of magni¬ 
tude change in the current for low biases. 

Besides an electronic switch, HMJs can serve 
as a thermoelectric mechanical switch, i.e. a 
junction who’s thermoelectric current - the cur¬ 
rent due to a temperature difference - can change 
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FIG. 5. Main panel: I— V curve of a 2,21-diaza[9]helicene (4) 
HMJ for two values of stretching distance, A r = — 2rA (solid 
blue line) and Ar = 2rA (dashed yellow line). A switching 
ratio of ~ 4 orders of magnitude between the ” On” and ” Off” 
state can be seen in the bottom-right inset, where the same 
is plotted on a log scale. Top-left inset: the thermo-current I 
as a function of temperature difference AT for the two states, 
showing a switching behavior. 

substantially. In the top-right inset of Fig. 5, the 
current is plotted as a function of temperature 
difference AT (at zero bias) for r = — 2rA (solid 
blue line) and r = 2 rA (dashed yellow line). We 
set the temperature of the right electrode to be 
room temperature, T rt = 300K, and the temper¬ 
ature of the left electrode is T L = 300 +AT. The 
thermo-current shows ~ 4 orders of magnitude 
difference between the two states at AT ~ 50K 
(in fact, since the thermopower changes sign, in 
principle one can have perfect thermo-electric 
switching because the ’’off” state can be tuned 
to have zero thermo-current). 

The dramatic change of conductance upon me¬ 
chanical stretching (or squeezing) is a direct pre¬ 
diction that is accessible to state-of-the-art ex¬ 
periments. In addition, since the whole effect 
is based on a coherent calculation of transport, 
it is of interest to examine the role of dephas¬ 
ing and decoherence in this system. Indeed, the 
majority of the theoretical approaches to trans¬ 
port and thermopower in MJs implicitly assume 
that the transport is coherent, and that there 
is no dissipation within the molecular junction 
itself, only at the electrodes or at the molecule- 
electrode interface [66, 67, 84], Recent experi¬ 
ments provide some evidence for coherent trans¬ 
port in the form of a zero-bias conductance dip 



A r [A] 


FIG. 6. Conductance as a function of the change Ar in the 
inter-electrode distance measured from the relaxed configu¬ 
ration, based on the DFT data (filled gray squares) and on 
the TB calculation, in the presence of dephasing, for different 
values of the dephasing time, Tdeph = 1CP 10 ,10 —11 ,..., 10 -1 's. 


[85-91]. However, these experiments typically 
compare between different junctions (measured 
in separate experiments), and lack a control pa¬ 
rameter that can be tuned to verify the presence 
of coherence in situ (except temperature in some 
cases [92]). As we demonstrate below, the de¬ 
pendence of conductance on electrode distance r 
can be an additional experiment to show coher¬ 
ent (or incoherent) transport in MJs. 

For this aim, we repeat the calculation of con¬ 
ductance vs. Ar (as in Fig. 4), now in the pres¬ 
ence of dephasing. We treat dephasing in a phe¬ 
nomenological way, by adding to the self-energy 
(Eq. 2) a diagonal term = —ir de ph, where 

r d eph = ft/ r deph and T dep h is the dephasing time, 
i.e. the average time it takes for an electron re¬ 
siding on the helicene molecule to loose its phase 
[67, 93-96]. In Fig. 6, the conductance as a func¬ 
tion of Ar is plotted for different values of the 
dephasing time, r dep h = 10 -10 ,10 _u ,..., 10~ 17 s 
(as in Fig. 4, the squares are the data from the 
DFT calculation, which is purely coherent). As 
seen, dephasing has a substantial effect on the 
conductance tunability, and dephasing times of 
~ 1CT 16 — 10~ 17 are in fact detrimental to the 
switching behavior. 



























FIG. 7. ZT as a function of Ar for a series of homolo¬ 
gous diaza[n]helicenes (from 2,15-diaza[6]helicene (1) to 2,31- 
diaza[14]helicene (9)). Inset: maximal thermoelectric FOM, 
ZT rnax . as a function of helicene length [n] (number of fused 
(hetero)aromatic rings). 


RESULTS: THERMOELECTRIC FIGURE OF 
MERIT 

With the TB parameters (and their depen¬ 
dence on inter-electrode distance r) determined, 
we can proceed to evaluate the thermoelectric 
FOM of HMJs. In Fig. 7, ZT is plotted as 
a function of Ar for diaza[n]hclicenes (l)-(9) 
differing in the number of all-ortho fused (het- 
ero)aromatic rings (n = 6, 7, ... 14) and, ac¬ 
cordingly, in their length. To calculate ZT real¬ 
istically, we take also the phonon thermal con¬ 
ductance into account, with a typical value of 
k = 20pW/K [10, 97-99]. Two conclusions can 
be quickly drawn from this figure: (i) ZT can be 
tuned by changing r, and displays order of mag¬ 
nitude difference within the same junction, and 
(ii) the longer the helicene the better its thermo¬ 
electric performance. We find that the optimal 
ZT is obtained at Ar ~ 0. In the inset the 
maximal ZT , ZT max is plotted as a function of 
helicene length (solid line), showing that ZT max 
values can be as large as ZT max « 0.6, which far 
exceeds current MJ values. 

To understand the origin of this strong 
ZT dependence on r we focus on the 2,27- 
diaza[12]hclicene (7) HMJ ([n] = 12) HMJ as an 
example (but the results are qualitatively similar 
for all [n]), and study how the different elements 
that constitute ZT depend on Ar. In Fig. 8(a) 


the transmission function T(E) is color-plotted 
as a function of E and Ar. The black arrow indi¬ 
cates the position of the LUMO peak, and it can 
be seen that upon compression of the junction 
(which increases the inter-stack hopping term 
1 1 ), the position of the LUMO changes, shift¬ 
ing down toward the Fermi level (dashed line). 
Once the LUMO crosses the Fermi, the HOMO- 
LUMO gap is no longer visible and the HMJ be¬ 
comes’’metallic” (resonant tunneling). The shift 
in the LUMO position implies that the filling of 
the gap cannot be viewed as due to the emer¬ 
gence of an addition conduction channel (the 
7r-stacking), but rather due to a change in the 
electronic properties of the helicene, induced by 
enhanced hopping along the stacks. 

The thermopower S, whose absolute value is 
plotted in Fig. 8(b), reflects the behavior of T(E) 
(since it is approximately equal to the logarith¬ 
mic derivative of T(E) at the Fermi level [4, 67]). 
The peak in |Sj at Ar ~ —0.5rA is a result 
of the Fermi level being close to (but not at) 
the LUMO resonance, where the slope of T(E) 
is maximal. Surprisingly, |»S'| increases dramat¬ 
ically when the junction is stretched, although 
the LUMO resonance shifts away from the Fermi 
level. This is due to the decrease in T, which 
typically increases the slope of the transmission 
function. Note that the thermopower can reach 
values of S ~ 100 — 200 pV/K (and larger for 
longer HMJs), which is more than one order of 
magnitude larger than values typically observed 
in molecular junctions [17-21], 

In Fig. 8(c) the ratio between the thermal 
and charge conductances n/oT is plotted as a 
function of Ar. In metallic systems, that ra¬ 
tio is expected to be ~ 1 in accordance with 
the Wiedemann-Franz (WF) law, and this is in¬ 
deed seen for the compressed HMJ, where the 
molecule is ’’metallic”. However, upon stretch¬ 
ing the molecule there is a significant deviation 
from the WF law. Thus, the dependence of ZT 
on Ar, plotted in Fig. 8(d), can be summarized 
as follows: the peak near Ar ~ OrA is due to 
the enhanced slope of the transmission function 
caused by the shift of the LUMO resonance to¬ 
ward the Fermi level. The peak at large Ar val¬ 
ues is due both to the increase in S (originating 
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Transmission function T(E) 



Ar [A] A r [A] 


FIG. 8. (a) Color-plot of the Transmission as a function of 
energy E and inter-electrode distance Ar. (b) Thermopower 
(absolute value) |Sj as a function of Ar. (c) The Lorentz 
number re/crT as a function of Ar. (d) ZT as a function of 
Ar. The calculation is done for 2,27-diaza[12]helicene (7) (n 
= 12, see text for other numerical parameters). 

in the decrease in T) and to the violation of the 
WF law. 


CONCLUSION 

In Summary, we have shown that hclicene- 
based molecular junctions can be used as a 
benchmark tool to investigate and mechanically 
tune transport and thermopower in molecular 
junctions. Using a combination of DFT and 
tight-binding calculations, we have shown that 
by compressing or stretching the molecular junc¬ 
tion, its electronic properties can be controlled, 
leading to drastic changes in the transport (from 
a so-called ”insulating” to ’’metallic” resonant¬ 
tunneling behavior) and possible substantial in¬ 
crease in the thermopower and thermoelectric ef¬ 
ficiency. The origin of these effects is the heli¬ 
cal spring-like structure of the hclicene molecule; 
due to this structure, stretching and compress¬ 
ing leads to changes in the tunneling matrix 
elements between vertically-neighboring atoms 
(Fig. 1). Control over the hclicene length and 
number of rings was shown to lead to more than 
an order of magnitude increase in the thermo¬ 
electric figure of merit of HMJs, and to a typ¬ 
ical thermopower which is an order of magni¬ 
tude larger than that of more common molecular 


junctions. The idea of a mechanically control¬ 
lable conductance and thermopower in molecular 
junctions should apply for other molecules with 
non-planar geometry, for instance cyclopara- 
phenylenes [100-102], cyclacenes [103], ball-like 
molecules [104, 105], carbon cages [106-108], tai¬ 
lored Fullerens [109], Fullcrene cages [110-112], 
and short DNA molecules (already demonstrated 
as elements in molecular junctions [113, 114], 
paving the way towards in-situ mechanical con¬ 
trol of transport properties of molecular junc¬ 
tions. 
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